Ice-lens formation and geometrical supercooling in soils and other colloidal materials. 
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We present a new, physically-intuitive model of ice-lens formation and growth during the freezing 
of soils and other dense, particulate suspensions. Motivated by experimental evidence, we consider 
the growth of an ice-filled crack in a freezing soil. At low temperatures, ice in the crack exerts large 
pressures on the crack walls that will eventually cause the crack to split open. We show that the 
crack will then propagate across the soil to form a new lens. The process is controlled by two factors: 
the cohesion of the soil, and the geometrical supercooling of the water in the soil; a new concept 
introduced to measure the energy available to form a new ice lens. When the supercooling exceeds 
a critical amount (proportional to the cohesive strength of the soil) a new ice lens forms. This 
condition for ice-lens formation and growth does not appeal to any ad hoc, empirical assumptions, 
and explains how periodic ice lenses can form with or without the presence of a frozen fringe. 
The proposed mechanism is in good agreement with experiments, in particular explaining ice-lens 
pattern formation, and surges in heave rate associated with the growth of new lenses. Importantly 
for systems with no frozen fringe, ice-lens formation and frost heave can be predicted given only 
the unfrozen properties of the soil. We use our theory to estimate ice-lens growth temperatures 
obtaining quantitative agreement with the limited experimental data that is currently available. 
Finally we suggest experiments that might be performed in order to verify this theory in more 
detail. The theory is generalizable to complex natural-soil scenarios, and should therefore be useful 
in the prediction of macroscopic frost heave rates. 



I. INTRODUCTION 

Frost heave is an incredibly powerful process, capa- 
ble of molding landscapes anywhere that temperatures 
drop below 0°C. During freezing, changes in soil volume 
cause the surface to swell and this can result in a vari- 
ety of fascinating patterns being sculpted into periglacial 
landscapes [H, Hf. These take many forms, such as stone 
circles, polygons and labyrinthine shapes 0, and occur 
on Earth as well as on the frozen surface of Mars 
Frost heave also leaves its mark after frozen soils have 
melted. Damaged rocks can break away to kickstart 
landslides, while topsoils can be washed downslope by 
meltwater. Even flat, frozen soils often collapse unevenly 
when melting, resulting in hummocked ground and dis- 
tinctive 'drunken forests' where trees fall, or are left at 
unnatural angles due to soil movement. This 'thaw weak- 
ening' is also well-known as a major cause of structural 
damage to roads and buildings in cold climates and has 
a significant cost, in particular because of the regular re- 
quirement to repair roads after freezing. For instance, 
in 1999 it was estimated that in the USA over two bil- 
lion dollars was spent on reparations to roads caused by 
frost-heave damage alone [5j- 
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In order to model the macroscopic effects of frost heave 
it is important to understand the micro- and meso-scale 
physics that underlie this process. However, despite ex- 
tensive prior research, there are still important questions 
remaining about how frost heave occurs. It is known that 
frost heave is not due to the expansion of water upon 
freezing (as commonly assumed); frost heave still occurs 
when water is replaced by substances that contract upon 
solidification [6]. Instead heave occurs when segregated 
ice forms within a column of soil that has its base in 
an unfrozen liquid reservoir. Thermomolecular pressure 
gradients then cause liquid to be drawn up from the rel- 
atively warm reservoir towards the cold, segregated ice 
whereupon it freezes (this process is commonly referred 
to as 'cryosuction') Q- Thus the soil 'expansion' is due 
to the added mass of water brought into the soil column 
from the reservoir (e.g. Q). 

In the last forty years there has been much experimen- 
tal work aiming to elucidate the mechanisms at work in 
the ice segregation process and it is now known that ice 
can segregate in a range of morphologies 0, [l(| ■ However 
the most well-known configuration is the periodic ice-lens 
formations that were first observed in the early experi- 
ments of Taber 0, [ll[ . It is the growth of these ice lenses 
that determines the soil's heave rate. Therefore it is im- 
portant for us to understand how they form in order to 
make accurate, quantitative predictions of the frost- heave 
process. Previous theories of ice-lens formation have been 
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proposed by O'Neill & Miller [ij], Gilpin pj] and Rem- 
pel et al. [8| which have provided the foundation for most 
subsequent work on frost heave (e.g. [H], El). However 
the theory of O'Neill & Miller is effectively empirical, as 
it is based on a postulated 'stress partition function' that 
is physically-unjustifiable. Rempel et al. Q provided a 
microphysically-correct model for the growth of periodic 
ice lenses. They, and Gilpin consider how periodic ice 
lenses can form in the presence of a 'frozen fringe' (a re- 
gion, ahead of the warmest ice lens, where liquid freezes 
into ice in the micron-sized pores of the soil) . However in 
some systems, such as colloidal clays, it is known that pe- 
riodic ice lenses can form when no frozen fringe is present 
(e.g. In order to explain these observations, we 

present here a new, physically-based theory that shows 
how multiple lenses can form, even in the absence of a 
frozen fringe. Moreover, we show that this model is in 
good qualitative and quantitative agreement with experi- 
ments on ice-lens formation in a variety of colloidal mate- 
rials. Throughout the paper, we shall refer to the mate- 
rial as a soil, however we note that this theory is equally 
applicable to general colloidal materials. 



II. A NEW MODEL FOR ICE-LENS 
FORMATION AND GROWTH 

In order to formulate a new model of ice-lens growth, 
we draw upon intuition gained from a series of experi- 
mental observations which indicate that the initial seg- 
regation of ice in freezing materials is seemingly 'crack- 
like' in nature. Firstly, Chamberlain & Gow [l8[ and 
Arenson et al. [l9| froze water-saturated soils and, by 
taking cross-sections of the frozen material, observed the 
presence of ice-filled shrinkage cracks in a polygonal pat- 
tern. They noted that these cracks are caused by the 
development of very low pressures (due to cryosuction) 
in much the same way that desiccation cracks form on 
the surface of rapidly-drying clays. Secondly, recent ex- 
periments have enabled us to closely observe growing ice 
lenses in order to investigate the nucleation process of 
subsequent lenses. Figure [T] shows three snapshots taken 
during the propagation of a new ice lens across a section 
of a freezing cell in our experiments on the freezing of 
kaolinite clay. The cell is cooled from above. The upper- 
half of each picture shows the last (warmest) ice lens, 
while the lower (light-colored) half is kaolinite clay. A 
new lens can be seen growing across the cell from right 
to left in a remarkably crack-like manner (the new lens 
grows approximately 1mm in front of the old lens, and 
the frames are taken 20 minutes apart). Similar observa- 
tions have been made in experiments on the freezing of 
silica microspheres by Watanabe et al. [HI (Watanabe, 
private communication). Thirdly, Akagawa et al. [20| 
performed multiple frost-heave experiments on a sample 
of Dotan silt. They concluded that the tensile strength 
of the soil strongly affects ice-lens formation. Ice-lens 
formation pushes soil particles apart, so that when the 



frozen soil is melted, there is a lower tensile strength at 
the points where previous lenses reformed. In subsequent 
freezing experiments, ice lenses formed at these weakened 
points. Interestingly, despite forming ice lenses at the 
same points, less undercooling was required for new ice 
lenses to form in later experiments. Akagawa et al. [2(| 
suggest that this is due to the fact that the forces acting 
to open an ice lens must overcome the tensile strength 
of a soil before an ice lens grows. Thus when the ten- 
sile strength is lower, it is easier to form ice lenses (we 
also note that ice-lens nucleation will be promoted by 
increased soil porosity that may exist at the position of 
previous lenses). Further evidence supporting this inter- 
pretation was obtained by performing freezing tests on 
the same material, broken up and reconsolidated in or- 
der to reduce the tensile strength of the soil. In this case 
inter-lens spacing was reduced by a factor of 5, suggesting 
that ice lenses are significantly easier to nucleate when 
the soil tensile strength is reduced. 

These experiments indicate the strength of the soil is 
an important factor, and that ice-lens nucleation may be 
caused by material fracture in the presence of internal 
stresses. Thus we take a fracture-mechanics approach to 
the problem. We consider the stresses on a small flaw in 
the material, and calculate at what point these stresses 
become sufficiently large to overcome interparticle cohe- 
sion and cause the flaw to enlarge. As we shall see, this 
will allow us to obtain a new, physically-intuitive condi- 
tion for ice-lens nucleation. We should note that fracture 
mechanics has been previously used by Walder and Hal- 
let [21( to model segregated-ice growth in freezing rocks, 
a process which has been recognized as being related to 
ice-lens formation in soils. However, we will take a differ- 
ent approach to that work, in particular by not assuming 
the presence of a frozen fringe. 

In order to clearly demonstrate the essential physics 
underlying our theory of ice-lens formation, we consider 
a simple, one-dimensional model of frost heave as shown 
in figure [5] It will be seen later that the model can 
be straightforwardly generalized to more complex set- 
tings. For continuity and transparency we base our 
notation on that of Rempel et al. The model 

is based upon laboratory frost-heave experiments us- 
ing directional-solidification cells such as have been per- 
formed by several groups (e.g. [l(| EH). A saturated soil 
or a colloidal suspension is placed in a cell in which a fixed 
temperature gradient VT is imposed. This 'frozen tem- 
perature' approximation is commonly invoked in model- 
ing direction-solidification experiments, and is based on 
having a cell with a large thermal mass in comparison 
to that of the material contained within it. Then tem- 
perature fluctuations caused by latent heat release on 
freezing, and due to the interphase variation in thermal 
conductivity have little impact on the overall tempera- 
ture gradient in the system (e.g. [Hj]). 

We define the vertical co-ordinate z so that z — when 
T = T m , the bulk melting temperature of pure water. 
The imposed temperature is below freezing for z > 0, 
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FIG. 1. (Color online) Ice lens formation in freezing kaolin- 
ite clay. The clay is being frozen from the top downwards in 
the apparatus used by Peppin et al. [22I ]. (i)-(iii) are pho- 
tographs taken 20 minutes apart in time-lapse footage of ice- 
lens growth. The dark, upper region of each photograph is 
the previous ice lens. The light, lower region is kaolinite clay. 
A new ice lens (dark, crack-like) can be seen propagating from 
right to left through the clay. A vertical shrinkage crack can 
also be seen on the right-hand side of the photographs. 
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FIG. 2. (Color online) Schematic diagram for the freezing of 
clay. 



and hence at some point above z = the soil will begin 
to freeze. In a typical experiment a sequence of ice lenses 
form across the cell. The position of the bottom of the 
warmest ice lens is given by zi, and the temperature at 
this point is XJ. We do not assume the existence of a 
frozen fringe and hence on the warm side of the warmest 
ice lens is a region of soil that may or may not be partially 
frozen. The reason that the soil can remain unfrozen be- 
low the melting point is due to the Gibbs-Thomson effect 
whereby the melting point of an ice surface convex into 
its melt is reduced relative to that of a planar ice-water 
interface [2j . This prevents ice from invading small radius 
pores until a certain undercooling of the soil is reached 
|8j. Regardless of the existence of a frozen fringe, pro- 



vided that there is an interconnected liquid network in 
the porous medium, we can measure the Darcy pressure p 
(see Appendix A for detailed discussion) . At some depth 
Zh < 0, we assume that there is a warm liquid reservoir 
at which the Darcy pressure can be assumed to be a con- 
stant, ph- In addition there is an overburden pressure P Q 
exerted on the upper surface of the soil (in field situa- 
tions this corresponds to the weight of material overlying 
the soil) and so mechanical equilibrium tells us that the 
ice pressure within any lens that extends across the cell 
is given by P a . 

A useful relationship between the Darcy pressure of 
pore water p and the pressure P of ice in local equilibrium 
with the water is given by the Clapeyron equation 



P -p = pL, 



{T m - T) 
T 



(1) 



Here p is the density of the solid/liquid phase (we ignore 
density differences between ice and water) and L m is the 
latent heat of melting of ice .24]. Care needs to be taken 
in applying this equation to ice in freezing porous media 
because strictly-speaking the equation is derived for liq- 
uid in equilibrium with bulk ice. In effect this means that 
it should be applied to ice that takes up a volume much 
larger than a particle size such as that in an ice lens or 
in a large flaw (see Appendix A for a discussion). Thus 
care must be taken in examining ice inside typical soil 
pores that are around the same size as a soil particle. As 
an example of the use of this equation we know that the 
pressure in an ice lens is P , as it is in mechanical equilib- 
rium with the overlying weight given by the overburden 
pressure. Therefore assuming that the liquid adjacent to 
an ice lens is in local bulk equilibrium with the ice lens, 
we find that the Darcy pressure of liquid at an ice lens is 
given by p = P —pL m (T m —T)/T m . We note that strictly 
the Clapeyron equation only holds when the growth rate 
of the ice is slow, and kinetic effects need to be added at 
higher growth rates. The implications of these effects for 
this model will be addressed in a forthcoming work. 

Here, we assume the existence of 'large flaws' within 
the soil in order to be able to use the Clapeyron equation 
((T|). By 'large' we shall take it that the length of the flaw 
is as least ten times larger than the radius of the con- 
stituent particles of the soil. Such flaws are comm only 
assumed to exist in the analysis of material fractures |25| , 
and they have been observed experimentally in a variety 
of clays. For instance Diamond [26| observed pore sizes 
of larger than 100/im in kaolinite, while scanning elec- 
tron micrographs showed individual clay particles to be 
O(l^m) in size. Similar results were obtained for ben- 
tonite, illite, and a range of silts. 

Within this framework let us consider the conditions 
under which a new ice lens will form. Consider a large 
flaw ahead of the warmest ice lens and at a temperature 
below T m , as shown in figure [3ja). We assume that ice is 
present within the flaw (we will return to the question of 
how this occurs in section lV Al) . As mentioned above, the 
flaw is assumed to be large compared to a soil particle of 
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size R so that we can make use of the Clapeyron equation 
(P), while we also assume that it is small relative to the 
distance to the nearby ice lens and to the size of the fully 
grown lens. The Clapeyron equation shows that, because 
the undercooling T m ~ T is positive, the pressure of the 
ice within the flaw, Pi, will be larger than the Darcy pres- 
sure p. Therefore as the undercooling becomes large, Pi 
will become large and positive, and so the ice within the 
flaw will exert a pressure on the walls of the flaw. Even- 
tually, this pressure will overcome the cohesive strength 
of the soil, and the flaw will extend. Later, we shall see 
that the flaw will grow parallel to the previous ice lens. 
Thus, as the flaw extends, it cuts across the width of 
the sample with the result being the formation of a new 
ice lens. This is shown schematically in figure E^a-c), 
with the vertical extent of the flaws being overexagger- 
ated for clarity (in practice the vertical opening will be 
much smaller than the length of the crack) . We note that 
there will actually be many such flaws, each at different 
orientations to the previous ice lens. However (as can 
be derived from the analysis of stress that follows later 
in this paper) the crack-opening stresses will be greatest 
on the flaws that are orientated parallel to the previous 
lens, and it is these flaws that will be the first to extend 
to form new ice lenses. Thus we only need consider flaws 
of this particular, 'horizontal' configuration. Also, note 
that we use a single, isotropic pressure to characterize 
the ice stresses in the flaw, rather than concerning our- 
selves with the possibility of differential stresses arising 
in the ice. In fact, as explained in Appendix A, because 
the ice nucleates from the liquid (which can maintain no 
differential stresses) we expect the ice in a flaw to remain 
under isotropic stress so that it can be characterized by 
the pressure Pj. 

We assume that the material around the flaw can be 
treated as a linear elastic solid. This has the benefit of al- 
lowing us to use linear-elastic fracture mechanics to inves- 
tigate the conditions under which the crack will extend. 
Unfrozen, overconsolidated soils of high clay concentra- 
tions have been observed to have substantial linear-elastic 
responses (e.g. [13 HI]) and so this is a reasonable first 
approximation to the material behavior. When there is 
pore ice in the soil (i.e. in a frozen fringe) its presence 
will further reduce viscous and plastic behavior of the soil 
and so the linear-elastic approximation is again expected 
to be reasonable. In general, the elastic moduli of the 
system will change with z due to changes in particle pack- 
ing fraction and ice fraction. However initial flaws will be 
much smaller than the lengthscale over which these mod- 
uli are expected to change. So in conducting an analysis 
of the stresses acting on the flaw, we can take the flaw as 
being surrounded by a homogeneous elastic material. 

As discussed above, the pressure Pi of the bulk ice in 
the flaw is related to the local Darcy pressure p by the 
Clapeyron equation (QJ. Therefore if we know the local 
temperature and Darcy pressure at the crack position, 
its internal pressure is P, = p + pL m (T m — T)/T m . In 
addition the far-field stress acting on the soil is given by 
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FIG. 3. (Color online) The opening of a flaw in a soil, (a-c) 
The process by which a large, ice-filled flaw in the soil ahead 
of a frozen fringe is pushed open to form a new ice lens, (d) 
Pressures acting on the system around the flaw, (e) Pressures 
acting on the flaw after linear superposition to remove the 
farfield stress. 



the overburden pressure P a . This means that the stresses 
acting on the system reduce to those shown in figure[3Jd) : 
an internal pressure acting to open the crack P^, and the 
overburden pressure acting to close the crack. We wish 
to ascertain whether this system of stresses can drive 
crack extension and thus ice-lens formation. Because the 
material is assumed to be linear-elastic, we can use linear 
superposition to show that this problem is equivalent to 
the case of a flaw under internal pressure Pi — P , as 
shown in figure [3je) [2o| . 

According to linear-elastic fracture mechanics the flaw 
will grow when the stress intensity factor Kj (the mag- 
nitude of the stress singularity at the flaw tip) exceeds 
the critical value Kj Cl a material property measuring the 
fracture toughness of the soil. In the case shown in figure 
Ete) the stress intensity factor is given by 



Kj = (Pi - P )yfa, 



(2) 



where a is half the length of the crack (e.g. [IH). There- 
fore the flaw will extend when 



P — P = 

1 7 ± n — 



Kic 
VttAP 



0-f 



(3) 



Here we have set a = XR, where A > 10 so that the 
size of the flaw is much bigger than the particle size. We 
also recognize that the right-hand side of this equation 
can be interpreted as the tensile strength of the soil <r t . 
Note that as the crack grows, the stress intensity factor 
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increases. This means that when the internal pressure 
is sufficiently large to cause a crack to grow, the crack 
growth will not halt. Instead as the crack grows, new 
ice will freeze in the newly-created crack volume, pres- 
surizing the walls of the crack with the ice pressure Pi. 
The stress intensity factor will then exceed Kj c (as the 
new crack is longer than the original crack) , giving rise to 
further extension. This positive feedback means that the 
crack will continue to grow through the material until it 
reaches the material boundary (catastrophic failure) . At 
that point, the ice- filled crack becomes the new warmest 
ice lens (figure Etc)). Additionally we may also note that 
the fact that larger cracks have larger stress intensity 
factors means that the largest cracks will break open to 
form ice lenses first. This means that in calculating the 
conditions under which a new ice lens will form, we only 
need to consider the growth of the largest flaws present 
in the system, and need not worry about smaller pores 
with A < 10. 

As with other fracturing events, crack growth is ex- 
pected to be rapid in comparison to the vertical growth 
rate of an ice lens. As an example, in our directional- 
solidification experiments on kaolinite clay, time-lapse 
footage suggests that the fracture extension velocity is 
at least a factor of 10 times faster than the ice growth 
rate. This factor is likely to be substantially larger but 
is currently constrained by the time between frames in 
the available experimental data. Therefore it is reason- 
able to pinpoint the position of the new ice lens to the 
point at which the stress intensity factor first becomes 
large enough that the crack starts to extend. This means 
that equation ([3]) will represent the point at which a new 
ice lens will form. This new condition can be compared 
to the previously-assumed condition that ice lenses will 
first grow when the interparticle pressure between soil 
particles first becomes zero or reaches a critical negative 
value (e.g. [IGllIl) 

and it will be seen that these are 
not equivalent. However, it is interesting to note that 
Akagawa et al. [2(| have previously suggested that a 
condition of the same form as equation © gives good 
agreement with their experiments on ice-lens nucleation. 



III. THE CRITERION FOR ICE-LENS 
NUCLEATION 



A. Geometrical supercooling 

Consider a parcel of liquid at temperature T and Darcy 
pressure p ahead of (warmer than) the warmest ice lens. 
Suppose a new ice lens were to be inserted next to the 
liquid such that the liquid and the ice are in local equi- 
librium. The pressure of the ice lens will be P a , and if 
the ice and liquid are to be in equilibrium they must be 
at the temperature Ti (the ice-lens growth temperature) 



that from equation ([1} is given by 



Ti = T m 1 - 



Po-P 

pL m 



(4) 



We can compare the actual temperature of the parcel of 
liquid T to the ice- lens growth temperature Tt. If the 
temperature is below the ice-lens growth temperature, 
then energy can be released if a new ice lens forms. Thus, 
in principle, in a soil with no cohesive strength a new ice 
lens can form wherever T <Tg. However, in cohesive soils 
ice lenses will not form spontaneously as soon as T = T% 
because formation of an ice lens is hindered by the diffi- 
culty of breaking apart the soil structure to accommodate 
the new lens. Hence it is possible for the temperature of 
the soil to be less than Tg without a new ice lens forming; 
this can be thought of as being supercooled relative to 
the ice-lens state. We shall call this geometrical super- 
cooling. We show in Appendix B that this is analogous 
to constitutional supercooling that occurs in the freezing 
of alloys and dilute particle suspensions. 



B. Nucleation criterion 

Now we return to the ice- lens nucleation criterion . 
Substituting equations ([!} (with P given by the pressure 
of ice in the flaw Pi) and ((4]) into ([3]) we obtain 



P i — P = 

1 l 1 o 



pL m (Ti - T) _ Kic 



Tn 



(•5) 



and with some rearrangement we find that the condition 
for ice-lens initiation is 



Tf-T 



T Kr 



TmCTt 



7 



(6) 



This reinforces the argument, outlined above, that in or- 
der for the soil to fracture and nucleate a new ice lens, 
there must be geometrical supercooling of the pore liquid. 
When this supercooling reaches a critical amount, given 
by the right hand side of equation ^ the soil breaks 
apart and a new ice lens is formed. Note that this criti- 
cal value is likely to depend on the initial packing of the 
soil particles (e.g. 1205), as this will alter the flaw sizes 
present in the material. As a check of intuition, it can 
be seen that when there is no cohesion of the soil around 
the crack, Kj c = and a new lens will form at the point 
when the pore liquid becomes supercooled. In fact, we 
shall see later that if Ti — T < and a new ice lens forms, 
the ice lens would melt, so in the absence of geometrical 
supercooling we would expect to see a single ice lens or 
needle ice at the freezing surface of a soil (primary heav- 
ing) Q. This suggests that the presence of geometrical 
supercooling is intimately linked to the appearance of 
new ice lenses in a freezing soil. 

Figure 3] shows the two representative situations that 
can potentially occur in a freezing soil. The figures show 
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FIG. 4. (Color online) Schematic diagram of the soil condi- 
tions that can occur under a growing ice lens. The continuous 
curves correspond to the ice-lens growth temperature, while 
dashed curves represent the actual temperature of the soil, 
(a) Geometrical supercooling occurs. Ice lens nucleation can 
occur at the first moment when the geometrical supercooling 
(AT max ) is sufficiently large for the fracture condition ([6]) to 
be satisfied. A new lens will form at this point (z„). This type 
of profile for Tg will typically arise during the diffusive evo- 
lution of the pore pressure field after rapid cooling at the ice 
lens, analogously to the appearance of constitutional super- 
cooling in alloys 29j. (b) No geometrical supercooling occurs 
below the ice lens and no new ice lenses can form. This will 
typically occur at slower rates of cooling. 



profiles of the temperature T and the freezing temper- 
ature Ti in the soil. For simplicity we assume that the 
temperature is a linear function of z, though this will 
not be true in general. At the bottom surface of the 
warmest ice lens [z = z{), the bulk ice of the lens is in lo- 
cal thermodynamic equilibrium with the pore water and 
so T = Tg. At the water reservoir position (z = Zh) the 
water is warmer than the bulk melting temperature T m 
and so T > Ti . This can be seen from equation §4§ along 
with assuming that the overburden pressure is greater 
than or equal to the reservoir pressure ph, as will gener- 
ally be the case in typical frost heave situations. Figure 
SJa) shows the case where \dT/dz\ > \dT(/dz\ and so ge- 
ometrical supercooling occurs whereas figure HJb) shows 
the case where \dT/dz\ < \dTg/dz\ and so no geometrical 
supercooling occurs. 

In the case where there is geometrical supercooling 
there will be a point at which the supercooling reaches a 
maximum value, Tg — T = AT max , as shown in^a). In 
general, AT max will increase as the freezing progresses. 
Therefore after some time AT max will exceed the critical 
supercooling T m Kj c / ' pL m \J nXR and a new ice lens will 
grow at the point of maximum supercooling, z ~ z n . This 
occurs a finite distance ahead of the previous lens, and 
so this provides a mechanism by which periodic, discreet 
ice-lens patterns can occur. 

New cracks propagate parallel to previous ice lenses 
because the internal pressure driving the crack growth 
is greatest at z = z n . Thus the greatest amount of en- 
ergy is released for cracks that continue to propagate at 



this height. In addition, such horizontal cracks have no 
shearing component, which means that the crack should 
propagate in a straight line without kinking out of plane 
[25j] . This agrees with the general behavior seen in typi- 
cal frost-heave experiments. Note that we only consider 
here an ideal one-dimensional system with no features 
such as shrinkage cracks or boundaries in the horizontal 
direction. It is likely that in the field and the laboratory, 
there will be more complex stress conditions in the soil 
due to soil imperfections, and interaction between the 
ice lenses and other features. We assume that these ad- 
ditional stresses will then help to explain the variations 
in ice-lens configurations that are sometimes observed ex- 
perimentally (e.g. (HI)- 

Qualitative conditions under which ice lenses form can 
be determined from the fact that for supercooling to oc- 
cur \dTg/dz\ must be large near the ice lens, and drop off 
as the soil becomes warmer (see figure H]). From equation 
O, P - p = P L m (T m - T)/T m and so 



dz 



T m dp 
pL m dz ' 



(7) 



while the flow of liquid in the soil is controlled by Darcy's 
law 



U 



k dp 
p dz ' 



(8) 



where U is the total volume flux of water through the 
soil, k is the permeability of the soil and p is the dynamic 
viscosity of water. Combining equations (0 and (jHJ gives 



dz 



pUT n 

pkL m 



(9) 



Considering the case where U is approximately constant 
throughout the depth of the soil column, it can then be 
seen that the only term on the right hand side of this 
equation that is not a constant is the soil permeability 
k. So when the permeability is small, \dTg/dz\ is large 
and vice versa. Therefore ice-lensing conditions are ex- 
pected to be found when k is small directly ahead of 
the warmest ice lens, and increases towards the water 
reservoir at z = Zh- There are three major situations 
in which this can occur. Firstly, if ice enters the pore 
space directly below the ice lens to create a frozen fringe, 
the permeability will be drastically decreased there due 
to the reduction in pore volume of water (c.f. Hansen- 
Goos & Wettlaufer [3(| , whose work also shows how even 
small amounts of impurities will play a large role in deter- 
mining the permeability after the formation of a frozen 
fringe). In particular, frozen fringes are expected to oc- 
cur in soils with larger pore sizes. Secondly, if the soil 
matrix is very compressible, particle rejection in front of 
the ice lens compresses the particles together leading to a 
reduction in porosity and permeability (e.g. [3l|). This is 
expected to occur in swelling clays like bentonite, which 
has been observed to remain unfrozen down to — 8°C [l^ |. 
Thirdly, desaturation of the soil adjacent to the growing 
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ice lens can occur due to the rejection of dissolved air at 
the growing ice lens and to the large suctions that can 
develop within the soil because of the low temperature 
(e.g. [l3,[13)- Desaturation was identified as an impor- 
tant factor for consideration in the frost-heaving system 
by Penner [ 33] , and its occurrence can lead to reductions 
in permeability by several orders of magnitude (3~ij . We 
have observed evidence of desaturation in our own exper- 
iments on the freezing of kaolinite clay, and will report 
on this in future work. Any of these three scenarios pro- 
vide a viable mechanism for geometrical supercooling to 
occur, and thus the occurrence of frost heave. 



IV. HEAVE RATE SURGING 

It is worth devoting some consideration to what oc- 
curs when a new ice lens nucleates, in order to make 
predictions of the change in heave rate during the pe- 
riod of formation of a new ice lens. Figures 03a) and 
(b) show the situation immediately before the growth of 
a lens, and immediately after the lens has grown across 
the width of the system. As soon as the ice lens grows, 
the geometrical supercooling vanishes at the position of 
the new ice lens so that the pore water is in local thermo- 
dynamic equilibrium with the ice in the new lens. This 
means that Tg — T at the new lens and the Tg(z) profile 
can be recalculated using equation @ , as we have shown 
qualitatively in figure [5f b) . From equation Q Tg is lin- 
early dependent on the Darcy pressure p, so the drop in 
Tg at the position of the new ice lens corresponds to a 
drop in the liquid pressure there. This means that the 
adjustment of the Tg profile can be thought of as the 
adjustment of the soil liquid to the new (lower) Darcy 
pressure at the ice lens. 

Equation ([9]) shows how the volume flux of liquid in 
the soil is related to the gradient of the liquid freezing 
temperature. Before the ice lens nucleates, the volume 
flux at the position of the new ice lens Uq = U(z n ) is 
continuous because dTg/dz is continuous (we assume k is 
a smooth function of z). However, as can be seen in fig- 
ure EJb), after the ice lens has grown across the system, 
dTg/dz is discontinuous and equation ([9]) shows that the 
volume flux is discontinuous with U(z+) = U\ < Uq and 
U{z~) = U2 > 0. This means that the flow upwards out 
of the lens is decreased while the flow upwards into the 
lens (from the reservoir) is increased. Therefore there is 
a net flow of water into the new ice lens that allows it 
to thicken. One implication of this is that there should 
be a surge in the heave rate of a column of soil (as mea- 
sured by the volume flux of water being sucked up from 
the reservoir) each time that a new ice lens is formed, 
followed by a gradual decay until the point that the new 
lens initiates. This was observed experimentally by Pen- 
ner [35[ and Akagawa et al. [2(| ■ It is interesting to note 
that the amplitude of the surge is controlled by the crit- 
ical supercooling T m Ki c j ' ph m \j 'ir\R, and hence by the 
fracture toughness of the soil. Thus it may be possible 




(a) (b) 



FIG. 5. (Color online) Schematic diagram showing the 
changes occurring in the soil upon ice nucleation. The dark, 
continuous curves represent the freezing temperature of the 
pore water, dashed curves represent the actual temperature 
of the soil, (a) Before an ice lens forms, AT max is sufficiently 
large that ice lens nucleation will occur at the point z n . Uo is 
the volume flow rate of liquid at z n . (b) After ice lens growth 
Tg becomes equal to T at the ice lens. This alters the flow 
rate at z n so there is a discontinuity between Ui and U2, the 
volume flow rates of liquid above and below the new ice lens 
respectively. As discussed in the text, the volume flow rates 
are proportional to the gradient of Tg. The light, continuous 
curve corresponds to Tg before lens formation occurs. 



to derive information about soil material parameters by 
measurement of the surge in heave rate. Additionally, no 
surge is predicted in a soil lacking cohesive strength. 

This argument also demonstrates the link between geo- 
metrical supercooling and ice lens formation, namely the 
fact that the formation of a new ice lens should not be 
possible in a freezing soil if there is no geometrical su- 
percooling of the pore water. If a new ice lens nucleates 
in the soil in which the pore liquid is not supercooled 
(e.g. figurelUb)), then an amended version ofOJb), tak- 
ing account of a non-supercooled liquid, will show that 
U\ > Uq and U2 < Uq so there is a net flow of liquid away 
from the new lens. This means that the lens will reduce 
in thickness and disappear. Note that although a new ice 
lens cannot form without geometrical supercooling, it is 
still possible that ice will grow into the pores of the soil, 
forming a frozen fringe. 

It is important to recall that we have assumed here that 
the soil is a rigid matrix so that this pressure adjustment 
is instantaneous, however for compressible materials the 
pressure adjustment will spread out diffusively from the 
position of the new ice lens over a finite time. The argu- 
ments presented above can be modified to take account 
of the transient pore pressure fields, and it will be seen 
that the qualitative conclusions are the same. 



V. QUANTIFICATION OF LENS GROWTH 
TEMPERATURES 

In order to demonstrate that this theory is quantita- 
tively consistent with observations, in this section we 
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estimate the temperature at which a new ice lens will 
form. We apply this estimate to two different materials 
for which experimental data is available: firstly, the silica 
microparticles used by Watanabe and coworkers [TH, |36[ 
and secondly, kaolinite clay, used in our own experiments 
on ice lens formation. In both of these cases, there exists 
experimental evidence that ice lenses can form without 
the presence of a frozen fringe [16l |17| . Thus pore ice is 
not expected to be present and so we use unfrozen values 
for the material fracture toughnesses. 

As described above, geometrical supercooling can be 
facilitated by reducing the permeability of the soil di- 
rectly ahead of the warmest ice lens. We have sug- 
gested three possible mechanisms: (1) The appearance 
of a frozen fringe. (2) Having a highly compressible soil 
matrix. (3) Desaturation of the soil ahead of the ice lens. 
Each of these three mechanisms has the potential to re- 
duce the soil permeability by at least an order of magni- 
tude. For instance, Burt and Williams [13] showed reduc- 
tions in soil permeability of over four orders of magnitude 
as a frozen fringe formed. Peppin et al. |ilOj summarise 
evidence that the permeability of bentonite changes by 
over four orders of magnitude as it compresses. Finally, 
Peroni et al. [38j report reductions in the permeability 
of kaolinite of over an order of magnitude upon desatu- 
ration. This means that the soil ahead of the warmest 
ice lens can be divided into two parts, as shown in figure 
[6j Directly below the ice lens, there is an affected region 
with very low permeability, while further away from the 
lens, there is an unaffected region with high permeabil- 
ity. Provided that the low permeability region is not too 
thin, we can assume that the majority of the pressure 
drop between the reservoir z = Zh and the ice lens z = z\ 
occurs in the low permeability region, giving a pressure 
profile like that shown in figure [5] (the pressure is ap- 
proximately constant in the unaffected soil region). This 
can be converted into Tg using equation Q , and plotted 
against T to determine the supercooling. Figure [B] shows 
that the maximum supercooling is expected to occur at 
the boundary between the affected and unaffected region 
of soil. At this point, p ~ Ph and so, making use of 
equation (@} the geometrical supercooling is 



T t -T = T„ 



1 



Pq -Ph 



T. 



(10) 



In the experimental systems described above, the free 
surfaces are open to the air so that P = ph = P a , where 
P a is atmospheric pressure. Thus 

Tt-Tn T m - T, (11) 

so the ice-lens nucleation criterion ([B]) becomes 

T m Ki c T m a t 



T 



■n 



(12) 



Using this equation with figures for typical flaw sizes and 
material fracture toughness, we can estimate the under- 
cooling T m — Tut which a new lens will appear. It should 



Ice lens 




Unaffected soil Sll P eI '- ^ 
cooling * 



T, 



FIG. 6. (Color online) Schematic diagram showing the su- 
percooling that occurs when the permeability of the soil is 
reduced in a region directly ahead of the ice lens. The af- 
fected region has a lower permeability due to the presence of 
pore ice, compression of the soil, or desaturation. The Darcy 
pressure is shown on the left hand side. On the right hand 
side, actual temperature (dashed curve) and ice-lens growth 
temperature (continuous curve) are shown. The maximum 
supercooling is approximately at the interface between unaf- 
fected and affected regions of soil. Compare with figure [4] 



TABLE I. Data specific to the experiments of Peppin et al. 
(3,CIj and Watanabe [3(|. Approximate fracture toughnesses 
for silica microspheres and kaolinite clay are taken from Allain 
& Limat [39l ] and Ayad et al. respectively. Particle sizes 
for silica microspheres and kaolinite clay are taken from [3(| 
and [4ll ] respectively. 



Constant 


Value 


Units 




273 


K 


p 


1 x 10 3 


kgm -3 




3.3 x 10 5 


J kg- 1 


Kic (Silica microspheres) 


1 x 10 2 


Pam 1 / 2 


Kj c (Kaolinite) 


1.5 x 10 3 


Pam 1 / 2 


A 


10 




R (Silica microspheres) 


l.lxlO" B 


111 


R (Kaolinite) 


lxicr 6 


111 



be noted that this will serve as a lower bound for the un- 
dercooling, as with a less extreme pore pressure distribu- 
tion than that shown in figure |6j the system will need to 
be at a colder temperature before sufficient geometrical 
supercooling is present to form a new ice lens. 

If we substitute values into this equation from table H] 
we find that in the silica microsphere system of Watanabe 
[3o| , new ice lenses are expected to form at an undercool- 
ing of T m — T 0.014°C. This figure is within measure- 
ment precision of the minimum undercooling he recorded 
of ~ 0.02°C (that is, the warmest temperature at which 
ice lenses were observed) . For kaolinite clay we find that 
new ice lenses are expected to form at an undercooling of 
T m — T w 0.22°C. Given the uncertainties in Ki c , A and 
R and the fact that our estimate represents a lower bound 
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for the undercooling, the theory is in reasonable agree- 
ment with experiments, in which lenses were observed to 
form at temperatures of approximately — 0.3°C. 

We are currently unaware of further data regarding ice- 
lens formation temperatures (for which the tensile prop- 
erties of the material are available) against which we can 
test our theory, and therefore we must look to future ex- 
perimental results to provide more rigorous validation. 
A simple verification would involve a systematic study 
of ice-lens formation temperatures in a variety of mate- 
rials with different tensile strengths. At a more complex 
level, particularly useful experiments would be those that 
were able to tune the tensile strength of a soil while leav- 
ing other parameters, such as permeability, constant. As 
an example this might be achieved by using a material 
consisting of colloids with tunable interactions, or by al- 
tering the salt concentration in typical soils (while being 
aware of other potential changes that salt will make to 
the system, e.g. [13] )■ Finally, it would obviously be 
of great benefit to be able to experimentally view the 
ice-lens formation process on the microscopic level, such 
as is routinely done for such systems in equilibrium 
and we hope that methods of achieving this will become 
available in the near future. 



A. The presence of ice ahead of the warmest ice 
lens 

A key part of our theory is the presence of ice in at 
least one large flaw ahead of the warmest existing ice 
lens, but there still remains the question of how ice can 
nucleate in suitable flaws. When there is a frozen fringe 
present, then ice can grow down through the soil pores 
to fill a flaw. However when there is no frozen fringe, 
there is no pore ice, and an alternative mechanism is re- 
quired. In this case, there are two ways in which ice can 
grow into a flaw: homogeneous nucleation (spontaneous 
ice growth), and heterogeneous nucleation (growth from 
a pre-existing source of ice or from the wall of the flaw). 
We can rule out homogeneous nucleation in typical clays, 
as the temperatures at which ice lenses form (typically 
~ to — 0.5°C, e.g. dH) are much warmer than the tem- 
peratures at which ice spontaneously nucleates in clays 
(— 6°C or colder [HI). Therefore it seems likely that ice 
must be nucleated heterogeneously, and we suggest here 
two manners in which this may occur. Firstly, if the 
soil particles are good nucleators of ice (such as kaolin- 
ite), then ice can be nucleated from the sides of the flaw. 
In this case it is likely that ice will first nucleate in the 
largest flaws in the soil, for the same reasons that large 
droplets of water containing heterogeneous particles will 
solidify earlier than small droplets (e.g. [43]). This will 
result in the situation shown in figure [3] 

The second method by which we suggest that hetero- 
geneous nucleation can occur is from the side of ice-filled 
shrinkage cracks that form ahead of a growing ice lens, 
as shown schematically in figure EJa) . These shrinkage 



cracks are caused by suction in the soil adjacent to the 
growing ice lens, which puts the soil into tension, caus- 
ing it to fracture in a manner analogous to the forma- 
tion of mud cracks in drying soils [181 ] . Similar to the 
case in drying soils, the shrinkage cracks in freezing soils 
generally take the form of a polygonal (in plan view) ar- 
ray of cracks propagating downwards through the soil, 
and these are commonly observed in tandem with ice- 
lens formation (e.g. @, [l8|, Gjl ) • These shrinkage cracks 
provide a ready supply of ice ahead of the warmest ice 
lens; if a suitably large flaw exists in the soil on the side 
of a shrinkage crack, ice can grow into the flaw to form 
an edge crack, as shown schematically in figure [TJb). A 
new ice lens can then nucleate through extension of the 
edge crack by the mechanism described previously. In- 
terestingly, we have seen evidence for this behavior in 
our directional-solidification experiments on the freezing 
of kaolinite, an example of which is shown in figure [1] 

When a new ice lens is nucleated from an edge crack, 
we need to make a small modification to account for the 
change of geometry from a center crack as shown in figure 
[3[ to an edge crack as shown in figure (TJb). The deforma- 
tion in the vicinity of a crack that breaks a free surface 
(an edge crack) is constrained less than that near a fully- 
embedded, center crack. Thus the crack opening for an 
edge crack is greater than that for a center crack for a 
given applied stress, and as K scales with the opening, a 
higher stress intensity factor results. A detailed analysis 
of this problem (e.g. [25j) then shows that stress intensity 
factor at the crack tip ^ is modified to become 

Kt = 1.12{Pi - P ) y/^E. (13) 

Since K for the edge crack is larger than that for a center 
crack, the critical condition for crack advance is reached 
at a lower value of internal pressure, which requires a 
lower degree of supercooling. This is seen in the modified 
supercooling condition, which becomes 

T £ -T= TrnKlc r (14) 

1.12pL m Vn\R 

Thus the critical supercooling for the edge crack is re- 
duced from the critical supercooling calculated earlier © 
by a factor of 0.89. This means that there will only be a 
small quantitative difference between the two cases. 

In the case of nucleation off the side of a shrinkage 
crack, we also need to consider how the presence of 
shrinkage cracks affects Tg through their effect on the 
Darcy pressure p. Shrinkage cracks appear in order to 
relieve the tensile stress that has built up in the soil, and 
then propagate steadily in front of the warmest ice lens. 
However, unlike ice lenses, which continue to thicken af- 
ter they have nucleated, the shrinkage cracks remain the 
same thickness (e.g. [19|]). This means that there is no 
liquid being transported to, or away from the sides of the 
shrinkage cracks, and so the horizontal pressure gradient 
can be taken as zero. Thus the pressure is approximately 
one-dimensional, and should not be significantly affected 
by the presence of shrinkage cracks. 
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FIG. 7. (Color online) Schematic diagram of ice-lens forma- 
tion in a compressible soil without a frozen fringe. Initially 
shrinkage cracks form and propagate downwards from the ice 
lens into the soil. New ice lenses can then nucleate from 
the shrinkage cracks in regions where sufficient geometrical 
supercooling exists. The solid and dashed lines show the ice- 
lens growth temperature and actual temperature respectively, 
with the new ice lens forming at the position of maximum su- 
percooling, as in figure [4] 



VI. CONCLUSIONS 

In this paper, we present a new, physically-intuitive 
theory of ice-lens nucleation. Motivated by experimen- 
tal observations of crack-like behavior during the growth 
of ice lenses, we have determined the conditions under 
which an ice-filled flaw will split open. Ice in the flaw 
exerts an outwards pressure on the walls of the flaw, and 
when this pressure exceeds a critical value, we have shown 
that the flaw will extend horizontally across the soil to 
form a new ice lens. 

We have introduced a new quantity, which we term the 
geometrical supercooling, that is intimately linked to the 
ice-lens growth process. We have shown that at any point 
in a soil, there is a certain temperature Tg (referred to as 
the ice-lens growth temperature) below which energy will 
be released if a new ice lens were to form. However the 
cohesive nature of the soil opposes its tearing to form a 
new lens, and so the soil can exist at temperatures below 
Tg without new lenses spontaneously forming. The geo- 
metrical supercooling, Tg — T, represents the cooling of 
the soil below this temperature Tg. We have shown that 
when the supercooling exceeds a critical amount, the ice 
pressure acting on the walls of an ice-filled flaw becomes 
sufficiently large to extend it, driving rapid horizontal 
growth to form a new ice lens. Furthermore, we have 
shown that the existence of geometrical supercooling is 
necessary for the formation of ice lenses. Our analysis 
suggests that if Tg — T < 0, although pore ice may form, 
a new lens cannot appear regardless of the ice-lens nucle- 
ation mechanism. This is because if a new lens were to 
form, pressure gradients would immediately be set up so 
that it would melt away. 

The geometrical supercooling necessary for the growth 
of ice lenses arises when the permeability of the soil di- 
rectly ahead of the newest ice lens is low in comparison to 
the permeability of the warmer soil away from the freez- 



ing front. We have suggested three main scenarios by 
which this may occur. Firstly, as is commonly assumed 
in frost heave theories, a frozen fringe may penetrate the 
pores of the soil directly ahead of the warmest ice lens. 
Then the ice in the pores substantially reduces the soil 
permeability to water flow. Secondly, if the soil is highly 
compressible (such as a swelling clay like bentonite) then 
particles will be compressed together ahead of the ice 
lens, reducing the pore space, and thus the permeabil- 
ity. Finally, desaturation of the soil can occur ahead of 
the growing ice lens. This may be either due to bubble 
formation caused by air rejection from the growing ice 
lens, or due to cavitation caused by the large suctions 
that develop in the liquid. The desaturated region will 
be much less permeable than the saturated soil further 
from the lens. The highlighting of these three scenarios 
by no means precludes the existence of other mechanisms 
by which geometrical supercooling will occur. However 
at least one of these three is expected to occur in most 
freezing experiments, which indicates that geometrical 
supercooling is relatively common in freezing soils. 

We have demonstrated that our theory is quantita- 
tively in line with the available experimental data by 
calculating the warmest temperature at which new ice 
lenses are predicted to appear in two separate systems, 
the freezing of silica beads and of kaolinite clay. We hope 
that focused experimental studies of the ice-lens growth 
process will provide data that can further help verify 
our theory, and to this end we have suggested a num- 
ber of experiments which would be helpful in analyzing 
the ice-growth process. Our theory is also in good quali- 
tative agreement with observations. The theory explains 
the crack-like appearance of ice lenses as they grow, and 
also explains the experimentally-observed surge in the 
heave rate as each new lens forms (e.g. [HJ). Impor- 
tantly, the theory explains the experimental finding that 
it is possible for discrete, periodic lenses to form even 
without a frozen fringe. We demonstrate that a partic- 
ular mechanism by which this can occur is linked to the 
growth of ice- filled shrinkage cracks into the soil ahead of 
the warmest ice lens. After shrinkage crack growth, ice 
lenses can nucleate off the side of the shrinkage cracks 
to form the type of patterns shown schematically in fig- 
ure [7] and observed in experiments on freezing clays (e.g. 
USUI)- Interestingly, in the case of ice lens growth 
without a frozen fringe, we have shown that the ice-lens 
formation can be predicted using only the unfrozen physi- 
cal properties of the soil such as permeability and fracture 
toughness, which are readily experimentally-measurable. 

In order to clearly demonstrate the physics at work 
we have restricted ourselves to a simple system with a 
linear temperature gradient. However the theory does 
not rely on this assumption, and can be readily general- 
ized to nonlinear temperature distributions in the soil. 
Furthermore, the condition for ice-lens nucleation can 
be straightforwardly extended to higher dimensions, and 
thus the model should be useful in analyzing more in- 
volved situations in order to predict heave rates and ice- 
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segregation. We hope that this work will be of interest to 
scientists across many disciplines, not only for research 
into the mesoscopic properties of freezing soil, but also 
due to its potential for up-scaling for the calculation of 
heave rates on macroscopic scales. For instance, we antic- 
ipate this to be useful in predicting the potential impacts 
of frost heave on engineering structures in cold climates, 
and also in exploring differential frost heave and the abil- 
ity of freezing soils to create complex patterns such as 
stone circles and patterned ground. Finally, it should be 
noted that although we have focused on the growth of 
ice lenses in freezing soils, this theory is applicable to the 
freezing of dense colloidal suspensions and other porous 
media. Thus this work may also be of use in the analysis 
of ice segregation in porous media such as gels and rocks, 
and in colloidal systems such as are used in the develop- 
ment of biomaterials and other materials with specialized 
engineering properties [iij ]. 

Appendix A: The Application of the Clapeyron 
equation to freezing soils 

In a porous medium, the Darcy pressure is defined 
as follows: at the point where we wish to measure the 
pressure, we connect the pore fluid to a large reservoir 
of water and control the pressure of the reservoir until 
there is no flow between the pore fluid and the reser- 
voir (i.e. they are in mechanical equilibrium). Then the 
Darcy pressure is given by the pressure in the reservoir, p 
(e.g. [43]). There are several benefits of using the Darcy 
pressure over the actual liquid pressure that would be di- 
rectly measured within the pores of the porous medium. 
Firstly, use of the Darcy pressure means that we do not 
have to worry about details of the geometry of the porous 
medium, or about the local effects of factors such as 
van der Waals and electrostatic interactions between the 
solid and liquid components in the pores. These factors 
may cause the local pore pressure to vary over sub-pore 
lcngthscales, making this pressure hard to define on a 
macroscopic scale. However the Darcy pressure is well- 
defined macroscopically. Additionally, the flow through 
the porous medium is related to the gradient in the pres- 
sure by Darcy's law 

k 

v = — Vp. (Al) 
A* 

Secondly, although first derived for flow through satu- 
rated porous media, Darcy's law still holds if air bubbles 
form or ice nucleates in the pores of the porous medium, 
provided that the permeability k is replaced by a new per- 
meability that depends on the connectivity and amount 
of the water remaining in the porous medium. This is 
often used in the study of unsaturated soils (e.g. [34] ]) 

Finally, by careful treatment of the thermodynamical 
problem the Darcy pressure within the soil can be re- 
lated to the ice pressure in certain cases. If we consider 
thermodynamical equilibrium between ice and water at 



the same temperature, the pressure difference between 
the ice pressure and the liquid pressure is given by the 
modified Clapeyron equation 

P l -p = pL ^-^ 1 (A2) 

m 

where we have ignored the small effects due to the density 
differences between ice and water (e.g. [13] )■ It is impor- 
tant to emphasize that this equation is derived for equi- 
librium between bulk phases. This is because for small 
volumes of ice or water, the effects of surface forces act- 
ing on the volume are important, and these will affect the 
equation above. For instance, consider a pore filled with 
ice in a liquid-saturated porous medium. The edge of the 
ice consists of ice/premelted film/solid interfaces where 
the ice meets particles (e.g. Q), and ice/water interfaces 
where the ice meets pore throats between the particles. A 
first effect is due to van der Waals and electrostatic inter- 
actions between the ice and the particles/water. There 
will be an influence of these forces on the pressure in 
the ice over a lengthscale S, which is typically small in 
comparison with the size of a colloidal particle 0. A 
second effect is caused by the liquid and the particles 
surrounding the ice each transmitting different stresses 
to the pore ice. This results in a non-uniform pressure 
distribution in the ice near the pore walls. For a suffi- 
ciently large pore, this pressure distribution will average 
out away from the pore walls leaving a uniform pressure 
field Pi in the middle of the ice, with the averaging oc- 
curring over a lengthscale comparable to the size of a soil 
particle R. Thus we speak of the bulk pressure Pi in an 
ice pore if the size of the ice-filled pore is much bigger 
than 5 and R. In effect, this means that ice in pores 
much larger than R in size can be thought of as bulk ice. 

Now consider an ice lens, such as that shown in figure 
[3J and the water in the pores directly adjacent to the lens. 
The lens is much larger than an individual particle size, 
and therefore can be treated as bulk ice. We connect up a 
reservoir to the pore fluid to measure the Darcy pressure 
and allow the system to come to local thermodynamic 
equilibrium. Then the bulk water in the reservoir is in 
equilibrium with the bulk water in the ice lens and so 
the conditions are satisfied for equation (|A2j) to hold. 
From mechanical equilibrium, we know that the ice-lens 
pressure is given by the overburden pressure. Thus 

P - P = P L^-I± (A3) 

J- m 

at the ice lens. 

Now as a second scenario, consider an ice-filled pore in 
a porous medium. Typically the size of the pore will be 
around the same size as the particles, and so we cannot 
treat the ice within the pore as bulk ice. However, if we 
have a large flaw that is much bigger than the individ- 
ual particle size, the ice inside the flaw can be considered 
as bulk (as mentioned in the text, there is evidence to 
suggest that these large flaws are present in most soils). 



12 



Then equation (| A2|) holds, and if we know the local tem- 
perature and Darcy pressure of the liquid adjacent to the 
pore, we can calculate the ice pressure Pi as with the case 
of the ice lens above. This is the basis of our calculation 
of the force acting to open a new ice lens. 

As a final note, we also need to justify our assumption 
that the bulk pressure in the ice is isotropic within the 
flaw (i.e., that the ice is not under differential stress). 
For a solid under differential stress, the melting point of 
a face of the solid varies depending on the stress at the 
face j!^]. A consequence of this is that if such a solid 
is placed into a liquid, it cannot be in thermodynamic 
equilibrium with the liquid. If one side of the solid is 
in equilibrium with the liquid, then another side under 
a different stress will have a different melting tempera- 
ture. This means that the differentially stressed solid is 
unstable with respect to the isotropically stressed solid 
and so any liquid and solid in equilibrium must both be 
in an isotropic stress sta te ( i.e., the stress can be char- 
acterized by a pressure) [46]. Because for slow freezing 
rates typical of frost-heaving soils pore ice is in equilib- 
rium with surrounding pore liquid which is characterized 
by the isotropic stress (Darcy pressure) p, the pore ice 
must be under isotropic stress, unless external stresses 
are later applied to the soil column. 



Appendix B: The relationship between geometrical 
supercooling and constitutional supercooling 

As some readers may be familiar with the concept of 
constitutional supercooling in the freezing of suspensions 
and alloys, it is worth demonstrating the links between 
constitutional and geometrical supercooling. Constitu- 
tional supercooling is a phenomenon that has been well- 
studied for its role in causing morphological instabilities 
in the freezing of alloys and suspensions of colloidal par- 
ticles [H, In a suspension of particles, the freezing 
temperature Tf of the suspension is given by 

P L m (r "'~ r/) =n(0), (bi) 

where II is its osmotic pressure, which depends on the 
volume fraction of solid particles <f> [9(. This equation 
implies that the freezing temperature is a function of par- 
ticle volume fraction only: Tf = Tf(cj>). If we know the 
temperature and concentration of a particle suspension, 
then the temperature T can then be compared to the 
freezing temperature, and if T < Tf, the suspension is 
said to be constitutionally supercooled. The name 'con- 
stitutional supercooling' therefore can be seen to come 
from the fact that the amount of supercooling is depen- 
dent on the constitution (i.e. concentration) of the parti- 
cle suspension (or solution in the case of freezing alloys) . 
A typical instance where this is known to occur is during 
the rapid freezing of colloidal suspensions. Colloidal par- 
ticles are rejected from the advancing ice interface and 



build up in a boundary layer ahead of the ice, and if 
freezing is rapid enough constitutional supercooling will 
occur. This manifests itself as a morphological instability 
of the ice/suspension interface (e.g. (Toj). 

The analogy with geometrical supercooling comes from 
the use of the definition that n = P — p for a colloidal 
suspension, where P is the overall, bulk pressure of the 
suspension and p is the Darcy pressure. If we insert 
this expression into equation (|B1[) . and recognize that 
the bulk pressure of a suspension is the equivalent of the 
overburden pressure P Q in a soil, we obtain 

By comparison with equation it can then be seen that 
Te is the generalization of Tf to porous media. Thus ge- 
ometrical supercooling is the extension of constitutional 
supercooling to porous media. 

Although there is the strong analogy between the two 
supercoolings, geometrical supercooling is a much more 
general concept than constitutional supercooling when 
considering colloidal systems. Constitutional supercool- 
ing is restricted to the case where a freezing suspension 
has a well defined osmotic pressure as a function of <j>, 
so that equation (IB1[) can be used. This means that it 
only applies to compressible suspensions of colloidal par- 
ticles, and cannot be extended to incompressible, packed 
soils. Furthermore constitutional supercooling only as- 
sumes the presence of one 'impurity' in the solute (in 
this case colloidal particles). On the other hand, geo- 
metrical supercooling is applicable to compressible and 
incompressible soils, and we have shown that it takes into 
account the presence of further impurities in the system 
such as air or pore ice. That pore ice can exist at the same 
time as supercooling is the main reason that we have cho- 
sen to call this phenomenon 'geometrical supercooling'. 
With constitutional supercooling, the suspension is su- 
percooled relative to the temperature at which ice can 
first appear, and the nucleation of ice will immediately 
cause ice to grow rapidly into the supercooled region, 
thereby removing the supercooling However, as we 
have discussed in this paper, in a porous medium with 
a cohesive strength, pore ice can exist stably while the 
system is supercooled. Thus it is important to note that 
the system is not supercooled relative to the temperature 
at which ice can first grow, as microscopic, pore ice can 
already be present. Instead the system is supercooled 
relative to the macroscopic ice configuration. 
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